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1. Introduction 

For more than 50 years our understanding of interacting electronic systems has been 
essentially based on the Fermi liquid theory. In this way of viewing electrons, the electron 
is essentially described by a (nearly free) plane wave excitation, with a well defined 
momentum and energy (see e.g. [1]). Most of the phenomena in solids can thus be 
understood in this wave description, with the crucial addition of the Pauli principle which 
embeds the antisymmetric nature of the total wavefunction of the system. The Pauli 
principle blocks states too deep in energy and ensures the existence of a Fermi surface. 
Such an approach has been extremely successful to explain quite remarkable phenomena 
such as superconductivity. It has been crucial also to interpret the subtle effects of 
disorder in solids. Indeed although the naive scattering of electrons on impurities can be 
understood in a semi-classical approach, more subtle effects like Anderson localization 
rests on the quantum nature of the particles and can be interpreted as wave interferences. 
Using both scaling theories and sophisticated field theoretical techniques, it is now known 
that free electrons are localized by disorder in one and two dimensions, whereas a mobility 
edge exists in three dimensions. [2, 3, 4, 5, 6, 7] 

However when the interactions become strong this "wave" starting point can be a 
poor description of the system. Indeed as was described a long time ago by Wigner 
[8], if the Coulomb interactions dominate the wave functions of each particle becomes 
localized in space, instead of being a plane wave. The electrons form a crystal (Wigner 
crystal (WC)). In such a crystal phase our intuition based on the wave description of 
electrons fails. The wavefunction of the electrons is localized so the electrons become 
discernable particles by their position, and the Pauli principle plays a much smaller role. 
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Quantum fluctuations are still there however since the vibration modes of such crystal 
are still quantized. One is thus facing a novel phase of quantum systems with quite 
fascinating properties to explore. For example, for the pure system if the interactions 
are reduced the quantum fluctuations will increase and one can expect the crystal to 
melt. What is the nature of this melting and to which phase it melts is of course a 
crucial, and still an open question. Whether this is simply a Fermi liquid in which the 
particle have recovered their "plane wave" nature or whether some intermediate phase 
exists is yet to be understood. Needless to say one can also expect quite interesting 
effects of the disorder in such systems. In the crystal phase the disorder effects cannot 
be viewed as scattering of waves but rather as the pinning of this crystal on the random 
potential due to impurities. It is thus interesting to know how this pictures of a pinned 
crystal connects with the standard description of interacting disordered systems, adding 
interactions on the top of disordered free electrons [9, 10, 11]. From a theoretical point of 
view this is quite a challenging question. Most of the theoretical approaches used for free 
electrons either fail or become much more complicated when interactions are included 
which makes it more difficult to obtain unambiguous answers. Perturbative calculations 
or renormalization group calculations can be made for weak interactions. Unfortunately 
they scale to strong coupling, which leaves the question of the large scale/low energy 
physics still open. The crystal description thus provides and alternative way to study this 
quite complicated problem. Even in the case of the crystal the competition of repulsion 
between particles, that wants the crystal to be nicely ordered, and the impurities that 
want the particles of the crystal at random position is a quite difficult problem. As a 
results, physics has been mostly computed in the limit of very strong disorder where 
particles are pinned individually [12]. A summary of the situation is presented in Fig. 1. 

In addition to this pure theoretical interest, there is also a direct experimental rele- 
vance. The progress of nanotechnology have made it possible to produce such quantum 
crystals. There are various situations where one can reach Wigncr crystallization. In 
order to enhance the interactions effects compared to the kinetic energy it is necessary 
to reduce the electron density (kinetic energy is killed faster than interactions at low 
densities). Thus semiconductors where one can reach very low density of carriers have 
been prime candidates to observe Wigner crystallization. Two dimensional electron gas 
(2DEG) where the density can be controlled by an external gate are good candidates. 
Because the crystallization requires quite low densities, one can help it by killing the 
kinetic energy of the electrons by application of a very strong magnetic field. In that 
case the electrons describe cyclotron orbits and their kinetic energy is quenched helping 
the interactions. This allows to reach Wigner crystallization for higher densities than in 
the absence of a field. Thus the first signs of Wigner crystallization were observed in 
systems where the fractional quantum hall effect was seen, but at higher magnetic fields. 

Some experimental data in these systems is shown in Fig. 2. This data prompts for 
an immediate question: in all these systems there is no direct evidence of the crystal 
structure. Some attempts have been made to image such systems but one still cannot 
resolve individual crystal sites [16]. This is at variance with classical crystals such as the 
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Fig. 1. - (top) combined effects of disorder D and interactions (denoted loosely U) in d = 3. 
Whether an intermediate phase exists between the crystal phase and a standard Fermi liquid 
phase is an open question, (bot) same thing in two dimensions. In d = 2 it is unclear whether 
a metallic phase exists, or whether a true crystalline phase can exist. Nevertheless, for strong 
interactions a description in terms of a crystal is a much better starting point than a plane wave 
description of the electrons. Although on these pictures the phase diagram is written in term 
of the interactions for simplicity, usually the density is varied or the kinetic energy is killed by 
application of a strong magnetic field in order to reach the crystal phase. 
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Fig. 2. - (bottom) Soundwave absorption by a two dimensional electron gas (2DEG) under strong 
magnetic field. The frequencies at which the sound is absorbed correspond to the eigenmodes 
of the crystal (see Sec. 2'2), and are interpreted as evidence of a Wigner crystal in this system 
(from [13]); (top Right) Transport properties of a 2DEG. At strong magnetic field an insulating 
phase appears, again suggestive of the formation of a pinned Wigner crystal (from [14]). (top 
left) Current vs. voltage characteristics. One clearly sees a threshold field needed to have 
conduction. This is again reminiscent of what one expects of a pinned crystal (Prom [15]). 



vortex lattice where good imaging techniques exists (see e.g. [17]), and thus where it is 
possible to directly "see" the crystalline order. Here one has really to infer the existence 
of the crystal from circumstantial evidence, i.e. from transport measurements (mostly) 
or some other indirect measurements (acoustic wave absorption, compressibility etc.). In 
order to know whether the transport experiments can be considered as a proof or not 
of the existence of the Wigner crystal is it thus specially important to have a reliable 
theory that allows to compute the transport properties. Such a task is far from being 
trivial given the complexity of the problem. 
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Fig. 3. - Metal-insulator "transition" observed at zero magnetic field in 2DEG systems. Is the 
insulating phase (low density) a pinned Wigner crystal ? (From [19]). 



Quite recently systems with even lower densities could be achieved. In these systems 
a "metal-insulator" transition was observed. There is considerable debate whether the 
metallic phase really exists [18], but in my sense a more interesting question is whether 
or not the insulating phase (which is the low density phase) is a pinned Wigner crystal. 
Resistivity data is shown in Fig. 3. Finally and although I will not insist on these systems 
in these notes, Luttinger liquids can be viewed also as an electronic crystal. A detailed 
review of the properties of such systems and their connection with disordered electronic 
crystals can be found in [20] . 

There is thus considerable interest in understanding the properties of these electronic 
crystals, which is the subject I will present in these lectures. It is clear that I cannot 
expect to cover such a subject in these few pages, so I will just give the basic ideas 
and some references. From the theoretical side, the concepts and methods used are 
quite different from the methods used usually to tackle disordered electrons, and I will 
introduce some of them. I will discuss the minimal model needed to describe an electronic 
glass in Sec. 2. I will present some basic concepts for disordered crystals in Sec. 3 and 
discuss also some incorrect preconceived ideas on the physics of such systems. In order 
to give a more precise solution I will present two more sophisticated theoretical tools 
in Sec. 4: the replica method and a functional renormalization group. Some of the 
results for equilibrium transport will be examined in Sec. 5. Finally some conclusions 
and perspectives can be found in Sec. 6. 

In complement and connection with the present notes there are three other recent 
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set of notes that present other aspects of disordered elastic systems and where further 
references to other reviews or original papers can be found. The notes [17] deal with 
classical crystals such as vortex lattices. [20] is a rather complete review on quantmn 
problems. It also presents other techniques and discusses in details one dimensional 
systems. The present notes are complementary with this review (there is some overlap) 
in the sense that I will much more insist here on two dimensional electronic systems. 
[21] deals with the specific case of Wigncr crystal, and discusses in much more detail the 
issues of compressibility and out of equilibrium dynamics that I will only briefly mention 
in these notes. 

2. — Basic Description 

2'1. Elastic description. - Starting from the full electronic Hamiltonian (fermions 
with interactions and disorder) is a near impossible task. In the crystal phase some 
simplifications can be made since the particles are now discernable by their position. This 
allows for a minimal phenomenological model to describe such a crystal [22, 23, 24, 21]: 
one assumes that the particles are characterized by an equilibrium position and a 
displacement Uj relative to this equilibrium position. In order to define uniquely the 
displacement one should not have topological defects such as dislocations in the crystal. 
I will come back to this point in Sec. 5'2. Prom the original quantum problem one has 
to define the "particles" of the crystal. If the wavcfunction is localized enough then one 
can indeed ignore the exchange between the various sites and thus define "particles" that 
have a size given by the extension of the wavefunction as indicated in Fig. 4. Of course 
the density fixes the lattice spacing a. These two Icngthscalcs (size of particle, lattice 
spacing) are independent and should be kept. One then has to deal with the vibration 
modes of this crystal. As usual one considers the interaction between particles 



and expands to second order in the displacements. If the potential V is short range this 
gives the standard elasticity. Here one has to be more careful since V is the (long range) 
Coulomb interaction between the electrons. One has thus to distinguish the two modes 
presented in Fig. 5. The compression mode changes the smeared density at one point 
of the crystal, creating an excess or a default charge. Such charges interact through the 
long range part of the Coulomb interaction. Since the density is essentially Vm (I will be 
more precise later), this mode gives a contribution V{q)q^u*Uci in the energy, where V{q) 
is the Fourier transform of the Coulomb interaction. Since V{r) ~ 1/r, V{q) ~ 1/q in 
two dimensions. The compression mode thus has an elasticity going as \q\u*Uq. On the 
other hand the shear mode does not change the density (averaged over a couple of unit 
cells), and thus the long range part of the interaction does not couple to this mode. The 
coupling comes from the short range part (or if one prefers V{q ^ 1/a))- The shear mode 
has thus still an elasticity going as q^ . Being more precise one can define the longitudinal 
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Fig. 4. - The three length characterizing the Wigner crystal. The size Ic of the "particles" in the 
crystal (at low temperature it is essentially given by the extension of the wavefunction around 
the equilibrium position, at large temperatures it is controlled by the thermal fluctuations and 
is the Lindemann length), a the lattice spacing is controlled by the density of particles, and the 
disorder is correlated over a length r/. The inset shows the triangular structure of the Wigner 
crystal. Particles are labeled by an equilibrium position and a displacement Ui.(From [24]) 



and transverse displacements for a mode q by 
(2) 



u{q) = -UL{q) + (- A zjuriq) 

q q 
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Fig. 5. - (left) Bulk compression mode. This mode changes the smeared density and is thus 
coupled to the long range nature of the interaction, (right) shear mode. This mode leaves the 
smeared density invariant and feels only the short range part of the interaction. 
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This leads to the action (see [24] for more details): 

(3) '^=9 / -^'^'^L{q,OJn){Pm(^l + CL{q))UL{-q,-OJn) 

+PmWcW„[ML(g,a;„)uT(-g, -Wn) - UT{q,UJn)uL{-q, -UJn)] 

denote the integration over the Brillouin zone Jg^ ^^^^2 and the cj„ are the standard 
Matsubara frequencies. The terms are coming from the /2m term in the Hamil- 
tonian, and are the ones embedding the quantum nature of the crystal. For a quantum 
crystal it is impossible to have the particles at rest since both its position and position 
would be determined simultaneously. Thus even at T = quantum fluctuations exists. 
The third term in (3) comes from the Lorentz force, that exists if the particles are sub- 
mitted to a magnetic field. The action corresponding to this force is simply the work of 
the force Fl ■ u where Fl = qdrU AB. Pm — is the mass density and Wc = e.B/m the 
cyclotron frequency. 

The precise coefficients CL,T{q) are the elastic coefficients for the longitudinal and 
transverse modes respectively. These coefficients can be obtained from an expansion of 
the coulomb correlation energy of the WC in terms of the displacements [25, 26]. As I 
discussed CL{q) ~ l^l and Criq) ~ q^- Keeping the lowest terms of the expansion one 
has 

CL{q) ^ d\q\ - bq^ 

(4) Criq) = cq^ 

where c,d,b are positive constants. 

2 '2. Consequences for pure system. - For a pure system the consequences of the 
quadratic Hamiltonian (3) are easy to carry out. The eigenmodes of the system are 
obtained by doing the analytic continuation ia;„ — > w + In the absence of magnetic 
field the longitudinal and transverse modes are decoupled. The longitudinal one has a 
dispersion 

(5) pmujl = CLiq) ~ d\q\ 

This is the standard plasmon dispersion in = 2. Because of the Coulomb repulsion it is 

extremely hard to excite this mode. The system is very rigid (much more than a normal 
elastic system) since it does not like to have its local density changed. This is not the 
case of the transverse mode which keeps a phonon-like dispersion 

(6) pmCOT = Criq) ~ cq^ 

In the presence of a large magnetic field the two modes are mixed. It is easy to obtain 
the new eigenmodes by diagonalizing the two by two matrix in (3). This gives the modes 
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zero field 



High magnetic field 



uj-{q) 



oc (trans.) 
oc q (long.) 



„3/2 



Table I. - g dependence of the eigenmodes in a Wigner crystal in the absence of magnetic field 
or for a very strong field. 



of Table I. The "plasmon" mode is now pushed to the cyclotron frequency lOc, while 
there is a low energy mode with a q^^^ dispersion. This low energy mode is the mode 
that was probed in the sound absorption experiment shown in Fig. 2. 

For the crystal, the current is simply given by J = ep^dtu, making thus the conduc- 
tivity very simple to compute since it is essentially the correlator of the displacements 
(up to a factor w^). The conductivity is thus given by 

(7) o-a/3(a;) = -{epof'i{ijj + i6){ua{q = 0,oJn)*Ui3{q = 0,Un))iu,r,^u,+iS 
In the absence of magnetic field the conductivity is simply 

(8) o-xx(w) - 



m u + iS 

The real part is thus just a ^-function at zero frequency. This traduces the fact that 
the crystal can slide freely in response to an external force (the electric field). The pure 
quantum crystal is thus a perfect conductor. In the presence of a finite magnetic field, 
the conductivity becomes 

/n\ r \ f \ ^^^0 -i{uj + iS) 

(9) (^xx{i^) = (^yy{l^) = 



m ~{uj + i5y- + cj^ 

Thus the (5-function peak is pushed at the finite frequency ujc- The d.c. conductivity is 
zero. This is due to the fact that in presence of a magnetic field the electrons describe 
cyclotron orbits, and do not drift in response to a static electric field. Paradoxically the 
system is then an insulator. On the other hand the Hall conductivity is given by 

(10) CF^y{Oj) - 



m —{u) + i6Y+uj1 
This gives for the Hall resistivity 

Oxv muc B 



(11) = 



which quite remarkably is the classical Hall resistivity. 
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Of course these results are for the pure system only and the crucial question, in order 
to make contact with the experiments, is to determine how the disorder changes the 
above results. 

2'3. Disorder. - The coupling of the crystal to the disorder writes 
(12) 5dis = JdrJ drV{r)p{r,T) 

The potential V represents the impurities. In these systems there are usually two main 
sources of disorder. One is provided by donors/acceptors far from the 2DEG. They 
interact with the 2DEG via the long range Coulomb interaction (see e.g. [27]). Since 
they are far from the 2DEG this produces a smooth potential. Another source of disorder 
is provided by interface roughness [28] . This type of disorder has a wavelength which is 
comparable or shorter than the lattice spacing of the crystal, and is most likely the main 
source of pinning. 

In the limit where the impurities are weak and the pinning is collective, the disorder 
can be modelled by a gaussian potential, (see [22, 17] for a more in depth discussion of 
this limit) 



(13) V{r)V{r') = Ar J {r-r') 

Aj-j, is a delta-like function of range r/ which is the characteristic correlation length of 
the disorder potential (see Figure 4). The density reads 

(14) p^r)=Y,S{r-R^-Ui) 

i 

where ^ is a 5- like function of range Ic (see Figure 4) and m = w(i?°). Since the disorder 
can vary at a lengthscale r/ a priori shorter or comparable to the lattice spacing a, the 
continuum limit Ui u{r), valid in the clastic limit \ui — ^ a should be taken 

with care in the disorder term. This can be done using the decomposition of the density 
in terms of its Fourier components [29, 30] 



(15) p(r) ~ po - PoV • u + po XI 



where po is the average density and K are the reciprocal lattice vectors of the perfect 
crystal. This is a decomposition of the density in Fourier harmonics determined by the 
periodicity of the underlying perfect crystal as shown on Fig. 6. Let mc examine the 
various terms. The first term is just the average density. The second one is the textbook 
expression of the density in terms of the displacements. It corresponds to a density 
averaged at a scale of a few lattice spacing (smeared density). This part of the density 
thus traduces local changes of the density due to compression modes similar to the ones 
in Fig. 5. The other terms describes change of the density at a lengthscale comparable or 
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Fig. 6. - Various harmonics of the density. If one is only interesting in variations of the density 
at lengthscales large compared to the lattice spacing one has the standard "elastic" expression of 
the density in terms of the displacements. In a crystal however on has to consider variations of 
density at lengthscales smaller than the lattice spacing and higher harmonics are needed [29, 30]. 



smaller than the lattice spacing. The sum over the reciprocal lattice vectors K obviously 
reproduces the 5 function peaks of the density (14). The finite range of 5 is recovered 
[30] by restricting the sum over K to momentum of order iiTmax ^ tt/^c- This expression 
of the density is quite useful since it allows to express the density at all scales in term 
of the slowly varying function u{r). The price to pay is the more complicated form of 
the operator that embeds in its very form the information about the periodicity of the 
crystal and the size of the particles. This is a quite useful trick, initially introduced 
for Luttinger liquids [31]. One keeps only the low energy- long wavelength part of the 
displacement field, since this field has a well defined continuum limit, and thus a quite 
simple continuum Hamiltonian (elastic Hamiltonian in this case). But it is still possible 
using this field that varies at a lengthscale large compared to the lattice spacing (this 
is the essence of the elastic limit) to describe phenomena that vary at a much shorter 
lengthscale, provided that one uses the proper form for the operator in terms of this 
continuum field. This expression of the density allows to write quite simply the coupling 
to disorder. Using the representation of density (15) it is easy to show [22, 20] that a 
Luttinger liquid is indeed very similar to the quantum crystals that are described here. 

Using the expression of the density (15) one sees that different Fourier component 
of the disorder play different roles. The Fourier components close to g ~ couple to 
Vu. These terms play little role (see [20] for more details). In our case because the 
compression modes are quite stiff thanks to the Coulomb repulsion this is even more so. 
The main effect of the disorder comes from the other Fourier modes, i.e. the ones that 
are close to q ^ K where K is one of the vectors of reciprocal lattice of the crystal. The 
effect of the disorder can thus be written as 
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Fig. 7. - The quantum crytal can also be viewed as a periodic array of lines that get pinned 
by columnar disorder, i.e. by disorder that is independent of one of the coordinates t which 
is the imaginary time for the quantum problem. The lines are the space-time trajectories of 
the particles. Although the action can be computed with the imaginary time r, to obtain the 
physical response for the quantum problem an analytic continuation is needed. This analytical 
continuation has an impact on many physical quantities, so the analogy between the classical 
system and the quantum one should be exerted with care. 



The coupling to the displacements is highly nonlinear which of course makes this problem 
very difficult to solve. The action (12) and (3) is the minimal model needed to describe 
a disordered quantum crystal. 

For a quantum problem, the disorder is time independent. Thus if one looks at the 
space-time trajectories of the particles, as shovirn in Fig. 7, one can view the system as a 
classical system of lines in d-f 1 dimensions [22, 20]. These lines are pinned by columnar 
defects, i.e. defects independent of one of the coordinates (r in this case). The problem 
of lines pinned by columnar defects has been studied extensively in the context of vortex 
lattices [32, 33], where in that case the defects are produced by heavy ions irradiation. 
Thus some features of the two systems are indeed similar and one can gain considerable 
insight by thinking in term of the classical system of lines. There is however an important 
difference. In the classical system t is simply one spatial direction, one is thus looking at 
the static of such systems. The dynamics of classical systems is a different problem, where 
an additional time t needs to be introduced (see e.g. [17, 34]). For the quantum problem, 
T is the imaginary time of the problem. Thus t dependence of the quantities is describing 
the equilibrium dynamics of the quantum problem. In order to get the proper quantum 
dynamics however one should make the analytical continuation r — > it (or more precisely 
iun — » Lu + iS). This analytical continuation, although seemingly innocent has severe 
consequences and changes many behaviors. Though although some analogy between the 
static (but z = T dependent quantities such as e.g. the tilt of the lines) behavior of 
classical crystals with columnar defects and the dynamics behavior of quantum crystals 
can be made, one should exert great care because of the analytic continuation needed for 
the quantum problems, which ultimately leads to some differences in physics. 
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3. — Basic description and preconceived ideas 

3'1. Relevance of disorder and basic lengthscales. - In order to define the basic prop- 
erties of such systems let me consider first what happens in the ground state of the 
problem. Since the disorder is time independent, one can consider all displacements as 
time independent to get the groimd state (it corresponds to the mode a;„ = 0). The 
problem is thus equivalent to a classical system in d dimensions. The time integration is 
just providing a factor (3 in front of the action and the problem is equivalent to determine 
the partition function of a systems 

(17) Z = j Vu[r]e-'^"'' 
where 

(18) = ^ / '^''^(^")' +P°j ^(^)P(^) 

I have forgotten for simplicity in the above formula the difference between longitudinal 
and transverse modes and simply rewritten everything with a standard elasticity (this 
corresponds in fact to the transverse mode with an elastic constant Cxiq) ~ 9^)- It is 
thus enough to examine a classical crystal, in presence of point like defects (for more 

details see [17, 34]). 

In order to understand the basic physics of such problem a simple (but ground break- 
ing !) scaling argument was put forward by Larkin [35]. To know whether the disorder 
is relevant and destroys the perfect crystalline order, let mc assume that there exists a 
characteristic lengthscalc Ra for which the relative displacements are of the order of the 
lattice spacing u{Ra) — u(0) ~ a. If the displacements vary of order a over the lengthscale 
Ra the cost in elastic energy from (18) is 

(19) iK-'a' 

by simple scaling analysis. Thus in the absence of disorder minimizing the energy would 
lead to i?a = oo and thus to a perfect crystal. In presence of the disorder the fact that 
displacements can adjust to take advantage of the pinning center on a volume of size 
iif allow to gain some energy. Since V is random the energy gained by adapting to the 
random potential is the square root of the potential over the volume i?^, thus one gains 
an energy from (18) 

(20) -VRi'^po 

Thus minimizing (19) plus (20) shows that below four dimensions the disorder is always 
relevant and leads to a finite lengthscale 



(21) Ra 
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at which the displacements are of order a. The conclusion is thus that even an arbitrarily 
weak disorder destroys the perfect positional order below four spatial dimensions, and 
thus no disordered crystal can exist for d < 4. This is an astonishing result, which has 
been rediscovered in other context (for charge density waves is known as Fukuyama- 
Lee [36] length and for random field Ising model this is the Imry-Ma length [37]). Of 
course it immediately prompt the question of what is the resulting phase of elastic system 
plus disorder? 

Since solving the full problem is tough another important step was made by Larkin 
[35, 38]. For small displacements he realized that (18) could be expanded in powers of u 
leading to the simpler disorder term 



where /(r) is some random force acting on the vortices. Because the coupling to disorder 

is now linear in the displacements the Larkin Hamiltonian is exactly solvable. Taking a 
local random force f{r)f(r') = AS{r — r') gives for the relative displacements correlation 
function and structure factor 



where () denotes the thermodynamics average and — denotes the disorder average. 
-Bthormai IS the samc correlation function in the absence of disorder due to thermal fluc- 
tuations (which remain bounded in d > 2 and are thus negligeable at large distance 
compared to the disorder term). Thus the solution of the Larkin model confirms the 
scaling analysis: (i) displacements do grow unboundedly and thus perfect positional of 
the crystal is lost; (ii) the lengthscale at which the displacements are of the order of the 
lattice is the similar to the one given by the scaling analysis. In addition the Larkin 
model tells us how fast the positional order is destroyed: the displacements grow as 
power law thus one can expect the positional order to be destroyed exponentially fast. 
However these two conclusions should be taken with a serious grain of salt. Indeed the 
Larkin model is an expansion in powers of u, and thus cannot be valid at large distance 
(since the displacements grow unboundedly the expansion has to break down at some 
lengthscale). What is this characteristic lengthscale? A naive expectation is that the 
Larkin model cease to be valid when the displacements are of order a i.e. at r = In 
fact this is too naive as was understood by Larkin and Ovchinikov. To understand why, 
in a transparent way, let me use the expression for the density (15) This immediately 
allows to reproduce the Larkin model by expanding the disorder term (18) 



(22) 




(23) 



B{r) = {[u{r) - w(0)]> = Sthermai(r) + -^r 



,i-d 



(24) 




K 



in powers of u. Clearly the expansion is valid as long as K^axU <^ 1 This will thus be 
valid up to a lengthscale Rc such that u{Rc) is of the order of the size of the particles Ic- 
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Note that this lengthscale is diflFerent (and quite generally smaller) than the lengthscale 
Ra at which the displacements are of the order of the lattice spacing a. The Larkin 
model cease to be valid way before the displacements become of the order of a and thus 
cannot be used to deduce the behavior of the positional order at large length scale. In 
addition it is easy to check that because the coupling to disorder is linear in the Larkin 
model, this model does not exhibit any pinning, since the model is invariant by changing 
u{r) u{r) + C where C is an arbitrary constant (remember that the total of the 
random forces has to be zero / drf{r) = 0). Thus sliding by an arbitrary amount does 
not cost any energy and any addition of an external force leads to a sliding of the vortex 
lattice. Thus the lengthscale at which this model breaks down is precisely the lengthscale 
at which pinning appears [38]. The lengthscale Rc is thus the lengthscale above which 
various chunks of the system are collectively pinned by the disorder. A simple scaling 
analysis on the energy gained when putting an external force 

(25) H = J d-^rF^^tuir) 

allows to determine the critical force needed to unpin the lattice. Assuming that the 
critical force needed to unpin the lattice is when the energy gained by moving due to the 
external force is equal to the balance of clastic energy and disorder one obtains 

for the critical force needed to move the lattice 

(26) i^c oc || 

This is the famous Larkin-Ovchinnikov relation which allows to relate a dynamical quan- 
tity (the critical force at T = needed to unpin the lattice) to purely static lengthscales, 
here the Larkin-Ovchinikov length at which the displacements are of the order of the 
size of the particle. Let me insist again that this lengthscale controlling pinning is quite 
different from the one Ra at which displacements are of the order of the lattice spacing 
at that controls the properties of the positional order. 

These two lengthscales Rc and Ra are thus the crucial ones in controlling the prop- 
erties of a pinned classical crystal, and thus will of course be crucial for our quantum 
crystals. For the quantum crystals this leaves us with two main questions: (i) what is 
the conductivity of the system, and which length controls it; (ii) what happens beyond 
the length Ra- 

3'2. Conventional wisdom and dynamics. - Up to know we have been on relatively 
firm ground, but answering the above questions, even for a classical crystal is a tough 
cookie. Given the difficulty in obtaining the answers, instead of a solution people have 
been happy for a long time with a kind of consensus on how the system "should" behave, 
based on physical intuition and some incomplete, or in some cases incorrect proofs. I 
will give in this section a summary of these preconceived ideas on disordered crystals, 
since they were at the root of most of the physical solutions given for the dynamics. 
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Fig. 8. - The (incorrect) physical image that was the commonly accepted view of what a dis- 
ordered elastic system would look like. The crystal would be broken into crystallites of size Ra 
by the disorder. Dislocations would be generated by the disorder at the same lengthscale. 



Fortunately recent progress have made it possible to truly analyze these systems. As 
usual the true behavior turned out to be quite different from what was naively believed. 
I will examine these solutions in the next section. But for the moment let me go along 
and describe the (incorrect) picture of a disordered crystal, that was used for a long time. 

The first point is to naturally assume that because there exists a lengthscale Ra for 
which displacements are of the order of the lattice spacing, positional order is destroyed 
exponentially fast beyond this lengthscale. This behavior would be in agreement with the 
too naive extrapolation of the Larkin model to arbitrary lengthscales and is in agreement 
with exact solutions of interface problems in random environments (so called random 
manifold problems) and solutions in one spatial dimension. It was thus quite naturally 
assumed that an algebraic growth of displacements was the correct physical solution of 
the problem, and thus that the positional order would be destroyed exponentially beyond 
the length Ra- This led to an image of the disordered vortex lattice that consisted of 
a crystal "broken" into crystallites of size Ra due to disorder. To reinforce this image 
(incorrect) "proofs" were given [39] to show that due to disorder dislocations would be 
generated at the lengthscale Ra (even at T = 0) further breaking the crystal apart and 
leaving no hope of keeping positional order beyond Ra- A summary of this (incorrect) 
physical image is shown on Fig. 8. 

If one assumes that all positional order is lost, this allows to solve for the dynamics of 
the problem. In that case one can consider that each crystallite is pinned independently. 
When one tries to move it it sees a pinning potential that one can approximate by an 
harmonic well. The shape of the well depends on the size of the crystallite. If the system 
is submitted to an electric field at frequency uj, each crystallite will respond as a particle 
in a well. Crudely speaking the minimal frequency one can excite is (using lo = cq and 

^min 

~ l/Ra) is 

(27) Wpin cx — 

The perfect conductivity of the pure crystal, i.e. the (5-function peak that existed at 



Electronic glasses 



17 




Fig. 9. - In an image where the system is broken into crystallites of sizer Ra, one can imagine 
that each crystallite is experiencing a pinning force and is in an energy well. There is thus a 
natural pinning frequency, ujp which is controlled by Ra- The conductivity has thus a peak 
at this frequency. Coupling between crystallites and distribution of pinning forces leads to a 
broadening of this pinning peak, very often chosen as lorentzian for lack of a better function. 



zero frequency is now pushed to the finite frequency Wpin, as shown in Fig. 9. If there 
is a finite magnetic field the peak at ujc is essentially unchanged and there is a second 
peak that appears at Wpin. Of course the approximation of independent crystallites is 
not exactly true and the coupling between the crystallites as well as the distribution 
of pinning strength leads to a broadening of the peaks. The broadening is hard to 
determine precisely since one does not know exactly the static solution, so quite often a 
lorentzian broadening is assumed. This physical image inspired from pioneering theories 
[36] used for charge density waves allows to compute, using some approximations, the 
optical conductivity (see e.g. [20] for further details on that point). This gives specific 
predictions, in particular the density dependence of the pinning frequency would be from 
the above mentioned approximations ujp oc -n}/"^ . 

Measuring the optical conductivity is thus a very powerful way to check for the pres- 
ence of a pinned Wigner crystal. An example of data is shown in Fig. 10. One clearly 
observes a pinning peak in the optical conductivity. However there are many problems. 
First the broadening is far from being simply lorentzian as is obvious when simply looking 
at the data. A much more serious problem comes from the density dependence of the pin- 
ning frequency that can be extracted from the above data and is shown in Fig. 11. When 
one compares the experiments, with these theoretical predictions the agreement is qual- 
itatively wrong. The predicted result would be totally opposite to the data which shows 
a decrease of the pinning frequency with the density. Such important problems when 
one tries to compare with the data could cast serious doubts on the interpretation of the 
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Fig. 10. - Optical conductivity for various densities for a 2DEG under strong magnetic field. 
The peak at a characteristic frequency (pinning frequency) is an expected characteristics of a 
pinned crystalline structure. The density and magnetic field dependence of the pinning peak 
provide a stringent test for the existence of a pinned Wigner crystal phase (from [40]). 
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Fig. 11. - Variation of the pinning frequency with the density for a 2DEG under a strong mag- 
netic field. The pinning frequency decreases with the density in all systems (different symbols), 
in contradiction with the naive calculations based on the physical image shown in Fig. 8. Such 
calculations would lead to uip oc v}^'^ in qualitative disagreement with experiments. (From [40]) 
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insulating phase in terms of a Wigner crystal and quite naturally other interpretations 
for this phase have been proposed [41]. 

In fact, the reason of this discrepancy is that the phenomcnological image of nearly 
independent crystallites is simply plain wrong. In order to be able to compute the 
transport properties we thus need sharper tools to solve the problem of pinned quantum 
crystals. These are the tools I now present. 

4. - Methods 

Let us see which tools of our theoretical toolbox one can use to tackle a disordered 
problem. If wc put a random term in our Hamiltonian. then every observable depends on 
the specific realization of the random potential. What is needed is then to take an average 
of these observables with respect of the probability distribution of the disorder to get the 
average response of the system. This is of course a theoretical trick. A real experimental 
system has usually only one realization of the disorder, and the self-averaging comes 
from the fact that the system is large enough so that little pieces of the system more 
or less behave as independent sub-systems. Solving the problem for a given specific 
realization of the disorder and averaging afterwards is of course a totally impossible 
program. In addition of being impossible it is in general stupid, since a given realization 
of the disorder breaks many symmetries of the system. Since a given realization of the 
disorder is not invariant by translation, all correlation functions depend now on both 
coordinates x and x' at which they are computed and not simply on the difference x — x' 
as in a translationally invariant system. On the other hand the averaged correlation 
function is invariant by translation so it is much simpler. Of course averages can be done 
order by order in a perturbation theory, but if one wants to go beyond perturbation one 
wants a method to average over disorder from the start. Unfortunately many techniques, 
that were useful for the non interacting systems (such as supersymmetry [7]) fail when 
interactions are included. I present here a quite general method that still works and 
is known as the replica method. The other useful method is the dynamical method 
(so called Keldysh technique) but it is of a more complex use. For a recent review on 
dynamical methods see [42]. This method is specially important if one wants to do out 
of equilibrium dynamics in quantum systems since this is the only one that can deal with 
this question. 

Let me assume that we want to compute the average value of some observable O for 
a system with a random potential V . The average value can be written as a functional 
integral over the configurations of the system as 



(28) 



{0)v 



/r»<^O(0)e-5^(<^) 



where 5^(0) is the action of the system for a given realization of the random potential 

V. Of course (0)y depends on V itself, so we have to average over all realizations of 
V. If we assume that the disorder has a probability distribution p{V) the average over 
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disorder is 

J'DVp{V){0)v 



(29) (O) 



JVVpiV) 



In general Sv {4>) is linear in V, something like (this is indeed the case in the disordered 
crystal problem) 

(30) Sv = So{(l>)+ J dxdTV{x)A{<t>{x,T)) 

Note that for a quantum problem the disorder is time independent. For the disorder 
one takes in general a gaussian disorder. This is very often justified by the central limit 
theorem. For example using a distribution 

(31) p{V) = e-^/'^^^^^)' = e'^^^^^'^" 
corresponds to the average 



(32) V{x)V{x') = D6{x - x') 

With the distribution (31) it would be very easy to perform the average (29) if it were 
not for the denominator in (28). Indeed in the absence of such denominator one has 

I / We"* / 'i^^(^)'e-/ dx J dTV{x)A{<j,{x,T)) ^ 

J'DVe~'^ I dxy(a:)2 J 



D 



(33) 



JdxJdTj dT'A{,l>{x,T))A{<t>{x,T')) 



One would end up with an effective action where the disorder has been eliminated and 
has given after average an interaction term so the action would be 

(34) 5eff = So{<f>) ^^JdxJdT JdT'A{cj>{x,T))A{cj>{x,T')) 

we could then treat this problem using our favorite method since it would not be more 
complicated than the type of problems that we already encountered in this book. 

The idea of the replica method [43] is thus to get rid of the denominator and to 
transform it into itor (sec the lecture notes by I. V. Lerner in this volume, for 

the use of replicas for calculation of the partition function). As with any really great 
method the idea is very simple. One can rewrite 



^^^^ JV4>e-SvW 



I 



V(f)e 



-1 n-l 

-SvW 
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with n = 0. If we forget n = for a moment and consider n as a positive integer 



n = 2,3,4,. 



then 



(36) 



-SvW 



I 



-Sv{<t>2 



-Svi4>n) 



where wc have introduced the fields (f)2, 4>z etc. The denominator can thus be rewritten 
as the product of n — 1 copies. The average (28) can thus be rewritten 



{0)v 



I>(?!>iO(^i)e-'^^('^i) 



(37) 



There is no denominator anymore. The price to pay is the introduction of n copies of the 
system. Of course one would only recover (28) if one can take the limit n ^ at the end. 
Before averaging over disorder in (37) all replicas (copies) are independent. Since there 
is no denominator in (37) one can do the average over disorder in the manner described 
above. One thus finds 

(38) ^ " / . . . P</'„0(</.i)e-^^« 
where the eflfective action is now 

(39) 5eff = X^5o(0a)-y J dx J dT J dT'AiMx,T))AiMx,T')) 



a=l,6=l' 



This is nearly the same form except that now one has n fields and the interaction couples 
all fields together. Wc have thus traded a disordered system with only one field for a clean 
interacting problem of n coupled fields. Of course the second one is more complicated 
due to the presence of the n fields, but as explained before we should be more equipped 
to tackle it. One important difficulty is of course that we should obtain a good enough 
solution for any n to be able to make the analytic continuation to n ^ at the end, 
since it is only in this limit that one recovers the disordered solution. Taking this limit is 
far from being obvious and contains hidden difficulties that I will briefly allude to later 
in this section. 

If we use this method we end up with an action S = S^i + S^is where S'ei is the elastic 
part of the action and the disorder part is 



(40) 



n p 

Sdis = -pQ^ X] X] / dxdrdr' cos{K{ua{r,T) - Ub{r,T'))) 



a,6=l K 
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The term (40) is of course extremely complicated to solve. To tackle such a term two 
paths are available. The first one rests on the hope that the properties of the system 
are dominated by small fluctuations aromid the minimum of one of the cosine. One 
would thus hope that expanding the cosine is relatively accurate. In fact things are more 
complicated because of the presence of replicas. The second possible way is to tackle 
the problem as a standard critical problem and try to find renormalization equations for 
the disorder term. Here again the situation is more complicated than one would naively 
think. I will examine these two methods in the following. 

4T. Variational method. — The first method is to try to reduce this problem to the 
"best" quadratic problem approximating 6*01 + Sdis- The most brutal way would be to 
simply expand the cosine, but as is well known for sine- Gordon problems this is very 
often a poor approximation. Here is would be a complete catastrophe for reasons that 
will become clear later. An improved way is to perform a self consistent calculation. The 
simplest way to implement such a scheme it through Feynman's variational approach 
[44]. The method has been extended to deal with replicas in [45]. Let me define a trial 
action So- Then the free energy F of the system verifies 



where ()o denotes an average taken with the action 5o. If So depends on some variational 
parameter one can hope that the "best" way (i.e. the one giving a physics as close as 
possible to the one of the true action S) of choosing this parameter is to minimize S'var- 
Let us take for So the most general quadratic action 



The choice of a quadratic action for So has two motivations: (i) quadratic actions are 
the only ones with which one is able to compute something; (ii) on a physical basis one 

can expect the properties of the system to be dominated by the small oscillations around 
one of the minima of the cosine. We will have to come back to this assumption. In 
doing so one implicitly neglects excitations that can take the system from one minimum 
of the cosine to the othc^r. For the sine-Gordon model these excitations are known as 
solitons [46]. For the pure problem this method is equivalent to the well known self 
consistent harmonic approximation, and gives excellent results. The disordered case is 
more involved [22, 20]. In particular special case should be paid to the mode a;„ = 
since the disorder is time independent. I will only sketch the main lines here. The full 
Green's functions are here the variational parameters. They are determined by 



(41) 



F<S, 



'var — 



Fo + {S- So)o 



(42) 




(43) 



'var 



= 



dGab{q,i^n) 



Electronic glasses 



23 



n n 



n 




Fig. 12. - (left) Replica symmetric structure. The matrix Gab has one value on the diagonal G 
and one value away from the diagonal G. (right) Replica symmetry breaking structure. The 
example of a two step RSB is shown. The matrix Gab is a hierarchical matrix. There is one 
value G on the diagonal and several values (here three) away from the diagonal. In the limit 
n ^ the value away from the diagonal becomes a function G{u) labelled by an index u £ [0, 1], 
giving its position in the hierarchical structure. In this example the function G{u) takes three 
values with two breakpoints (hence the name two step RSB). 



Using (41) leads to nasty integral equations, roughly of the form 

(44) G-,\q,cUr.) = {u;l + q')Sat + J rfr[l - cos(c.„r)]e-^ S.^^. . . . 

One has to remember that the limit n ^ should be taken at the end so Gat is in fact a 
0x0 matrix. In order to invert this matrix to solve the equations (44) one needs to know 
what is the structure of such matrices. In the original action, Sd is purely diagonal in 
replicas, whereas Sdis couples each replica with another and is thus independent of the 
value of the indices a, b. This suggests that a good structure for Gab is to take one (for 
each value of q and w„ of course) value on the diagonal, and one unique value away from 
the diagonal 

(45) Gab = [G - G]Sab + G 

Such a matrix structure is shown in Fig. 12. This is the natural choice, and is called 
replica symmetric since it respects the permutation symmetry between the replicas. The 
matrix depends in fact on only two numbers: the value on the diagonal G and the off 
diagonal value G. It is easy to see that such matrices can be inverted and have a well 
defined limit when n ^ Q. Calling G^^ and G^^ the diagonal and off diagonal value of 
the matrix G~i^ and the connected part Gc — Gab one has the inversion formulas (in 
the limit n ^ 0) 
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This allows to close the integral equations (44) and in principle to solve them. Unfortu- 
nately, it is easy to see that (45) cannot be right. This can be seen if one brutally expand 
the cosine. This indeed corresponds to a RS structure and gives "directly" So (this ap- 
proximation corresponds to the variational approach if one assumes the RS structure 
and that the disorder A is large). In that case only the a;„ = frequency appears due 
to the double integration on r and r'. Since the action is quadratic this means that the 
disordc;r would only affect the mode Ci;„ = 0, and thus would be unable to change the 
current. This is obviously completely opposite to the idea of pinning, so such an approx- 
imation misses all this physics. The answer to this paradox resides in the assumption 
of the replica symmetric structure for the matrix Gab- Indeed as is now well known for 
spin glasses, such a saddle point might be unstable [47]. In that case one should look 
for another structure for 0x0 matrices. This structure was discovered by Parisi, and 
breaks the replica symmetry [48]. The matrix is now a hierarchical matrix. There is 
not a single clement off diagonal but many values defined in a hierarchical way. I will 
not enter in the structure of the theory, and refer the reader to [49, 20] for a discussion. 
A pictorial representation of the matrix is shown in Fig. 12. Such matrices can also be 
inverted in the n limit. This structure is called replica symmetry breaking. In this 
limit the off diagonal clement becomes a function G{v) of a continuous variable varying 
between zero and one. If G{v) is a continuous function one speaks of continuous RSB, 
whereas if G{v) has step(s) one talks about n-step RSB (n being the number of steps). I 
will not dwell too much on the physical interpretation and refer the reader to [20] . As is 
well known for spin glasses, the occurrence of RSB signals glassiness, and the presence 
of metastable states. Since there is one matrix Gab{tl,^n) for each value of q and a;„, 
the replica structure depends a priori on q and utn- For quantum problems since the 
disorder is time independent one can show [22] that one has off diagonal replica terms 
only for the mode w„ = 0. The structure is replica symmetric for q > l/Rc and RSB 
for q < l/Rc- This confirms the physical picture that pinning and metastability appear 
above the Larkin length Rc. 

4"2. Functional renormalization group. - RSB might look like magic("'^), so let us try 
to solve our disordered problem by a more conventional renormalization approach. Doing 
it on quantum problems is tedious so I will explain the method on a classical system. 
One simply takes ?l = in the quantum action which removes all time dependence in 
time of u(r, r). The system is thus described by the classical action 



where I have retained only a single harmonic in the disorder term for simplicity. To 
implement an RG procedure one would decompose Ua{r) into fast and slow modes. If for 
example A is a cutoff in q space and A' a smaller cutoff, one defines (f2 is the volume of 



(47) 





(^) This is metaphysical question, since the only point is whether a method works or not! 
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the system) 



Wa(r-) = u<{r)+u>{r) 



(48) = ^ E "»(«)^""^ + ^ E "'^W^""^ 

q<A' A'<g<A 

One then integrates over the fast modes (A' < g < A) to get an effective action with new 
coupHng constants for the slow modes. Since the cutoff has changed it is then convenient 
to rescale space by g = q'A'/A or r = r'A/A', such that the new space variable q' has 
again a cutoff A. One can then iterate this procedure. Even without performing the 
integration over the fast modes let me look at the bare dimensions of the operators (tree 
level). Let us assume that u scales as u ^ _L'». One sees directly that periodic systems 
are peculiar. Because the disorder term is a cosine, with a periodicity which corresponds 
to the lattice spacing, this periodicity cannot be changed in the renormalization. This 
imposes C = [29, 30]. Let me parametrize the cutoff by A(Z) = Aqc"', where Aq is the 
initial cutoff. In order to keep the elastic part of the action invariant, one should thus 
have 

11/ A(l) " ''-^ 

(49) - ' ^ ' 



T{l + dl) Til) \A{l + dl) 

the d — 2 comes from the integration d'^r and the V^. This gives the scaling dimension 
for the temperature 

(50) ^=Til){2-d) 

We thus see that for d > 2 the temperature scales to zero. This is in agreement with the 
fact that the competition between disorder and elasticity is an energy problem, where 
the temperature plays little role. Let us look at the renormalization of the disorder term. 
Performing the same rescaling of d'^r and using (49) for the temperature gives 



A(; + dl) 



d-i 



(51) A(/ + d/)=A(Z) 
giving the equation 

(52) ^^ = A(0(4-rf) 

We thus see that the disorder is relevant below d = A. Note that around d = A the 
temperature is indeed irrelevant. One could thus naively think that this problem is 

quite similar to the standard 4 — e renormalization in critical phenomena. There is 
however an important difference. Let us compare the present problem with a standard 
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(non disordered problem) in the standard (j)"^ theory. Let me denote the action for this 
problem 

(53) ^ J d''r{V<f>r + d'^rcf,^ + d'^rct>^ 

In that case and contrarily to the disorder case both the elastic term and the "potential" 
terms have the same power of the temperature. This is of course the signature that a 
normal phase transition results from a competition energy-entropy. The extra factor of 
1 /T in front of the cosine in (47) is the way for the replicated action to indicate that before 
averaging one had to compare the two energies of the elastic term and of the disorder 
term. For (53) one absorbs the factor e~(''~^)''' coming from the spatial rescaling in the 
rescaling of u and one leaves the temperature invariant. This traduces the fact that a 
problem at T = Tg stays at Tg. The rescaling for u thus becomes 

(54) (j,'^{l + dl) = (l)'^{l)e^'^-^^'^^ 
This gives for the renormalization of the coupling constants 

(55) ^ = ^(6 -2d) 

Thus all terms in the potential that have high (i.e. larger than 4) powers of u are 
irrelevant. Thus regardless of the precise form of the potential, one can always expand in 
power series and the only important term is the cj)'^ term. The situation is totally different 
for the disordered system [50] . In that case u does not scale. It means that all terms in a 
power series expansion of the potential have the same scaling dimension and are equally 
relevant around four dimensions. Thus instead of having a single RG equation for a 
single parameter A one has now an infinity of (a priori coupled) RG equations for all the 
coefficients of the terms in the power series expansion of the potential. In other words 
for the disorder problem one has now to renormalize the whole function corresponding 
to the potential term. This means that one has to start with a potential term 

(56) Sais = E / ^'^^KW - Mr)) 

ab 

where A{z) is a function whose initial {I = 0) value is A;=o(-z) = cos(i4rz) but which will 
change form in the course of the renormalization. Because one has to keep the whole 
function under renormalization this procedure is known as functional renormalization 

group. I will not derive the equations to the next order (a good description of the 
derivation can be found in an appendix of [51]) but just give the result (for a classical 
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crystal) [29, 30] 



(57) ^ = eA{z) + lA"{zf - A"{z = 0)A"{z) 

Needless to say this is much more complicated to solve than if one had a single parameter. 
But in addition there is an important feature of the solution. Because we eliminate modes 

Ua becomes smoother and smoother thus in a sense Ua — ui, becomes smaller and smallcir 
under renormalization. Thus if the function A stayed analytic one could expand it in 
powers of u to get 

(58) A(n„ - Ub) = A(0) + ^A"(0)(«„ - ut,f + ^AW(0)(t.„ - ut,f + ... 

and one would be asymptotically dominated by the (ua — Ub)^ term. This would give back 
a RS solution, and in fact exactly Lar kin's result. Since we know that this is incorrect 
beyond a certain lengthscale, something has to go wrong with the above expansion. The 
only way to escape this is if A('')(0) becomes very large (in fact infinite) so that the 
above expansion breaks down and that higher power continue to play a role. Thus to 
get the glassy physics the function A(z) should become nonanalytic beyond a certain 
lengthscale. One can check that this is indeed the case. The equation for A^^' (0) is from 
(57) 

(59) = eAW(O) + 3[A(4)(0)]2 
which indeed explodes at the lengthscale 

/ 6 



(60) e' = 



V3Ac 



This corresponds indeed to the Larkin scale (compare with (21)). Beyond this lengthscale, 
although the function itself remains finite, it develop a non-analyticity. If one looks at 
the second derivative of A it develops a cusp, as shown on Fig. 13. I will stop there on 
the FRG method since this would take us way beyond the level of these notes. For more 
on the FRG and further references see [49, 34, 52, 53] 



5. — Bragg glass and disordered Wigner crystal 

51. Elastic model and experiments. - Let me now examine what are the results for 
the disordered Wigner crystal. I will not reproduce the calculations here that can be 
found in [22, 20] for quantum crystals in general and in [23, 24] for the specific case of 
the Wigner crystal. 

Using the variational method on (3) and (12) one finds that the most stable solution is 

an RSB one. This is suggestive that the pinned system indeed has glassy properties due 
to the pinning. I will thus use loosely the term Wigner glass in the following. Second, 
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Fig. 13. - In the FRG approach one has to renormalize the fuU correlator of the disorder. One 
starts with an analytic function (left). For crystals this function is essentially a cosine. After a 
certain lengthscale, corresponding to the Larkin length, the function develops a nonanaliticity 
(cusp) at the origin (right). This nonanaliticity is the signature in this method of the glassy 
properties of the system. 



since we know have from the displacement-displacement propagator, one can compute the 
conductivity. An example of the results is shown in Fig. 14. At first sight the agreement 
with the experimental data is much better (see [24] for a full discussion). Note that the 




Fig. 14. - An example of the conductivity of the Wigner glass as computed by the variational 
approach. Note the non lorentzian bradening of the peak and the magnetic field dependence of 
the pinning frequency (from [24]). 
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broadening is non lorentzian in agreement with the data. More importantly one finds 
a decrease of the pinning frequency with the density as ujp oc n~^/^, as well as a good 
magnetic field dependence of the pinning frequency This very good agreement gives a 
good confirmation that the insulating phase in the 2DEG under strong magnetic field is 
indeed a Wigner crystal collectively pinned by impurities. 

What is the reason of this difference with the naive answers explained in Sec. 3'2? 
An inspection of the variational solution [23, 24] shows that in the full solution the 
pinning frequency is not controlled by Ra as in the naive answer but by Rc (in fact the 
formula is more complicated). This is a very satisfactory and physical result. Indeed 
as we discussed, for a classical crystal Rc is the length that controls the pinning force. 
Although is it of course impossible to directly extrapolate it is quite reasonable to obtain 
that the pinning frequency, that is another manifestation of pinning, is also controlled 
by this lengthscale. This difference is crucial since, as in the Wigner crystal Ic and a are 
quite different Rc and Ra corresponds to quite different Icngthscales and have in general 
quite different dependence in the various parameters. Since Rc depends on the size of 
the particle Ic, this gives, for example for the case of strong magnetic field for which Ic is 
just the cyclotron orbit, additional magnetic field dependence to the pinning frequency. 
This is to be contrasted with charge density waves for which Ic ^ a due to the nearly 
sinusoidal density modulation and thus Rc ^ Ra- Thus borrowing directly the solutions 
that have been developed [36] for this case is dangerous and gives part of the physics 
incorrectly. 

I refer the reader to [23, 24], for other physical properties such as the Hall conductivity. 
I will just mention here the compressibility. Defining the compressibility for a charged 
system is a delicate question since one has to maintain the neutrality of the system. 
Naively one relates the compressibility to the density- density correlation function by 

(61) K = lim {p{q, Un = 0)p{-q, a;„ = 0)) 

Let me look first at the pure system. The compressibility is simply (only the longitudinal 
mode plays a role) 

(62) K{q) = lim 

If only short range interactions arc present in the system the longitudinal mode is a 
phonon-like mode CL{q) oc q^ and the one recovers a finite compressibility. On the other 
hand if one has long range Coulomb interactions Ci(q') oc \q\ and the compressibility 
becomes zero. This is simply due to the fact that (61) measures the density response to 
a change of chemical potential while keeping the neutralizing background unchanged. A 
charged system thus does not remain neutral, hence the infinite compressibility. One has 
thus to define the compressibility more precisely. Many derivations of the compressibility 
use instead directly a derivative of the free energy with respect to the number of particles. 
The free energy can be computed for a neutral system for an arbitrary number of particles 
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Fig. 15. - Capacitance measurement, which gives access to the compressibihty of the system. A 
voltage difference ^ is applied to a capacitor. Here for simphcity the capacitor is made of two 
plates of the 2DEG. 



which solves the above-mentioned problem. Unfortunately very often the calculation is 
only possible in some sort of approximate way such as an Hartree-Fock approximation. 
Here again the link with the direct measurements of the compressibility is not clear. 
Using such procedures, so called "negative" compressibilities are found for some range of 
the interactions, for interacting electrons. Similarly, experiments measure such negative 
"compressibilities" [54, 16] (see [18] for further references and discussion on this question). 

In order to make the physics of such negative compressibility more transparent, I 
discuss now a very simple way to compute it [21]. It is in fact a direct calculation 
of the quantity that is actually measured to determine the "compressibility" , i.e. the 
capacitance of a system made by with the 2DEG [54]. For simplicity I take here a 
capacitor formed of two identical systems, as shown in Fig. 15. Taking one system and 
one metallic plate would not change the results in an essential way. The Hamiltonian of 
the system is thus 

(63) f/?+i/2"+ E / Jnr-r')[p.(r)-po][p/3(r')-po] + ^ /[Pi(r)-P2(r)] 

If one assumes that the system is neutral in the absence of /i, then the charge on one 
plate when a potential jjL is applied is 

(64) (pi) = |[(pipi) - {piP2)] 

in linear response. (64) give directly the capacitance {pi)/fi. 

One can use the general formula (64) to compute the capacitance for the Wigner 
crystal. One substitutes in (3) the density decomposition (15). The Vu terms give 
directly the contribution of the long range part of the Coulomb interaction 

(65) -fflong-rangc ^PoXl [^"/3 (^X (q)^/? (-?)] 

q a/3=l,2 
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Since Vii{q) ~ 1/g, (65) gives obviously the part proportional to q in the elastic coeffi- 
cients for an isolated plane. The higher harmonics give the regular part (i.e. the part 
proportional to in the elastic coefficients. Such a way to determine the coefficient is 
equivalent the calculation of the coefficients in [25]. 

Taking a pure system the Hamiltonian becomes (only the a;„ = term of the longi- 
tudinal part needs to be computed to have the compressibility) 



{66)H 



4(9) \ ^ f cl^il) + pll'V,^iq) ph'VMq) \ ( u\{-q) 

,{q) J V pWVM cl^{q)+plq^V^,{q) J\ul{-q) 



where c|^(g) is the "short range" part of the elastic coefficients. Using (64) and the 
expression of the density for small q from (15) pL{q) = —poquhiq) one gets for the 
capacitance ^ = -^j^ h where Cgeom = l/(47rrf) is the standard geometrical capac- 
itance and 

(67) J- = hm ^ 

Cel 9^0 plq^ 

(again the factor of 2 comes from the fact that here I took two identical plates). The 
electronic capacitance corresponds to the propagator where only the short range part of 
the elastic coefficients is kept. Using [25] c^{q) = —ujQ{0.18..){aq)'^, where lvq = ^^^^3 , 
one finds for the Wigner crystal a "negative" compressibility. The fact that a system 
of discrete charges can lead to such effects has been noted before for classical Wigner 
crystals (see e.g. [55] and references therein). 

The disorder can be considered easily. The variational method discussed in Sec. 5 
gives straightforwardly the density-density correlation function which is simply related 
to the (mm) correlation. Since it is a gaussian approximation one can easily use the 
capacitance method shown above. The capacitance is thus given by (67) but where 
one should use the propagator in the presence of disorder. It is given by (within the 
variational approximation used [23, 24]) by: 

(68) (uu) ^ 



PmUjl + CL{q) + S(l - 6„fi) + I{0Jn) 

where S and I{cOn) are respectively a constant and a function related to the disorder 
verifying 7(0) = 0. At = (68) leads to a compressibility in identical to the one of 
the pure system, and thus also "negative". More on this point can be found in [21]. 

5"2. Defects. - We thus sec that our elastic description of an electronic glass is a 
very good starting point to extract the physical properties since we have now reliable 
analytic methods to tackle such problems. However, in view of the commonly accepted 
point of view that disorder would induce topological defects in a crystal (as explained 
in Sec. 3'2) is it important to ascertain the validity of the elastic description itself. In 
fact the creation of defects by disorder has been grossly overestimated and the elastic 
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Fig. 16. In d = 3 a classical crystal plus disorder is in a Bragg glass state (for weak disorder). 
Although the disorder destroys the perfect positional order of the crystal, it does it smoothly. 
The topological order of the crystal is preserved and there are no unpaired defects such as 
dislocations generated by the disorder. The system is thus nearly as ordered as a perfect crystal. 
In particular the relative displacements of two neighbors remain small, and thus such a system 
is totally described by an elastic theory. 



theory is in fact much more stable to the presence of topological defects than initially 
anticipated following the naive picture of Sec. 3'2. 

In d = 3 it is now known that below a certain threshold of disorder no topological 
defects can be induced by the disorder. The disordered elastic system is in a Bragg 
glass state [30] with a quasi long range positional order, much more ordered than the 
incorrect image of Fig. 8 suggests (see e.g. [17, 34] for a discussion and references on 
this point). The distance between unpaired defects is thus infinite in this case. Though 
although the crystal is distorted by the disorder and the perfect positional order is lost, 
the deformation is smooth and the topological structure of the crystal is preserved. The 
displacements of two neighbors remain small even if large distortions can accumulate at 
large lengthscales. The elastic description is thus perfectly accurate at all lengthscales 
in that case. An image of the Bragg glass is indicated in Fig. 16. 

In d = 2 the situation is marginal, and defects do appear in the ground state, but 
at distance much larger (for weak disorder) than the lengthscale Ra and not at that 
lengthscale [30, 56], as shown in Fig. 17. Since the ratio ^o/Ra can become arbitrarily 
large for weak disorder, it means that all important effects of disorder (they occur at the 
lengthscales Rc and Ra) can be reliably described by an clastic theory. Since the pinning 
peak is controlled by the length Rc <^ Ra ^ it means that dislocations will spoil 
the results of the elastic theory for the optical conductivity only well below the pinning 
peak, as shown on Fig. 18. This implies that the theory is a reliable tool to compute the 
characteristics of the peak and above, and thus most of the a.c. transport. Thus optical 
conductivity measurements are an accurate way of probing the Wigner glass nature of 
the system, since one has a reliable theory to compute the characteristics of the pinning 
peak in a Wigner glass phase. This allows for a detailed testing, as I showed above 
for the case of the 2DEG under strong magnetic field. For the systems in the absence 
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Fig. 17. - In d = 2 for a classical crystal, defects are generated by disorder at large lengthscales, 
so an elastic description is not valid beyond the distance between these topological defects. 
However for weak disorder is much larger than the lengthscale Ra (right). This is in marked 
contrast with the incorrect naive picture of Fig. 8 (left), where both length would be of the same 
order. It means that an elastic theory can reliably be used to describe the effects of disorder 
below (from [56]). 



of a magnetic field shown in Fig. 3, a measurement of the optical conductivity could 
thus provide an unambiguous answer on the presence of a Wigner glass in the insulating 
phase. Such an experiment remains to be performed. 

On the contrary this shows that to know d.c. conductivity it is necessary to take 
topological defects of the crystal into account (vacancies, interstitials, dislocations etc.). 



0)p 0)p 




Fig. 18. - (a) If dislocations occurred at scale Ra and the pinning frequency was controlled by 
Ra, as was naively believed, the elastic theory is incapable of giving any reliable information 
on the pinning peak. It would be necessary to include dislocations from the start, (b) As was 
discussed in the text, dislocations occur in fact at S> Ra and the pinning peak depends 
on Rc <C Ra- Thus the pinning peak is given quantitatively by a purely elastic theory. It 
is necessary to take into account topological defects such as dislocations only at much lower 
frequencies, and in particular if one wants to obtain reliable results for the d.c. transport (from 
[24]). 
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This makes the calculation of the d.c. conductivity much more difficult since at the 
moment there is no theory able to take into account topological defects, elasticity and 
disorder at the same time. It is thus impossible to make specific predictions specific 
to the Wigncr crystal and that could be quantitatively checked. The only case where 
such a calculation of d.c. conductivity can be done reliably is the case of one dimension 
(Luttinger liquids) [57]. 

6. — Further steps and unsolved questions 

In these short notes I have discussed some of the physics of quantum particles in 
the presence of strong repulsion. In that case the physics changes drastically from the 
standard "wave" description of indiscernable particles. The particles localize to form 
a quantum crystal. Examples of such systems are the Wigner crystal and Luttinger 
liquids. In such a quantum crystal, quantum effects are still important but particles 
become essentially discernable by their position, making the problem much simpler to 
tackle. The effects of disorder in such a disordered crystal are also quite different than 
on noninteracting electrons, and can be viewed as the pinning of the quantum crystal on 
impurities, in a similar way than for a classical crystal. 

I have briefly described the methods used to tackle such disordered elastic systems. 
Using them the a.c. transport and thermodynamic quantities such as the compressibility 
can be reliably computed. A pinned quantum crystal exhibit glassy properties, in a way 
similar to the Bragg glass for classical systems (quantum Bragg glass). There is good 
agreement with the predicted a.c. transport and the observed behavior of a 2DEG under 
a strong magnetic field, making a strong case for a Wigner crystal phase in such systems. 
The compressibility is found to be negative both for a pure Wigner crystal and in the 
presence of disorder, and detailed comparison with experiments on that point would be 
clearly fruitful. One important point is that the elastic description is an excellent starting 
point to obtain the a.c. properties of such systems. Indeed, although in rf = 2 topological 
defects in the crystal (dislocations etc.) are generated by disorder, the distance at which 
unpaired defect appear is very large compared to the scale at which the effects of disorder 
manifests themselves. This makes the a.c. transport a probe of choice to decide on the 
presence of a Wigner glass phase, since direct imaging techniques are not yet at the point 
where the crystalline structure can be seen. 

The d.c. transport is a much more difficult question because of the need to explicitly 
take the presence of topological defects into account. There are however some answers 
that can be obtained when the crystal is set in motion at finite velocity. In that case 
the defects become less important, and the properties are again crucially dependent on 
the periodicity of the system (Moving glass phase). As was shown in [58, 59], periodic 
moving periodic systems have a quite specific dynamics. Indeed due to the existence of 
the periodicity in the direction transverse to the direction of motion, the motion cannot 
average completely over the disorder. The moving system is thus submitted to a random 
potential, which leads to a channel like motion. The very existence of such channels has 
an important consequence if one tries to make the system move in the transverse direction. 



Electronic glasses 



35 



Indeed although the particles do move along the channels, the channels themselves are 
pinned. This means that even above the longitudinal threshold field, if one tries to 
apply a force in the transverse direction a transverse critical force still exists. For the 
two dimensional electronic system, putting a magnetic field is a simple way of applying 
a transverse force. If the lattice is sliding at velocity v, it is submitted to a Lorentz 
transverse force Fl = evB. The existence of the transverse critical force Ftr thus implies 
that the channel structure should not slide as long as < i^tr- There will thus be no 
hall voltage generated. On the other hand when Fl > -Ftr the channel structure should 
slide and a Hall voltage exists. The periodicity of the crystalline structure thus implies 
that one needs a finite longitudinal current before a Hall voltage exists. Experiments [60] 
have shown that a finite longitudinal current is clearly needed to develop a Hall voltage, 
in good agreement with the existence of a transverse threshold. The existence of such 
an effect is another direct probe of the crystalline (existence of a transverse periodicity) 
nature of the phase. 

We see that there are thus efficient, if not easy, ways via transport to check for 
the presence of a pinned crystalline phase (electronic glass). Most of the experiments 
discussed in the preceding sections were for a 2DEG under a strong magnetic field. It 
would of course be extremely interesting to use the same techniques to probe the 2DEG 
in the absence of the magnetic field, and analyze the experiments in the line of what was 
discussed above to check for the existence of a Wigner crystal in these systems. Among 
the interesting possible experiments one can note: 

• Measurement of the optical conductivity. In particular the density dependence of 
the pinning peak can be directly checked against the theoretical predictions of the 
pinned Wigner crystal. 

• I did not dwell too much of the non-linear transport and on the existence of a 
threshold field (as shown in Fig. 2), but a threshold field is also a characteristics of 
a pinned crystal. Since it is controlled by the Larkin length (26) it has to be related 

to the pinning frequency. So a comparison between the threshold field in the d.c. 
transport and the pinning frequency is another probe of the pinned Wigner crystal. 

• The Hall tension vs the longitudinal current (i.e. the measure of the transverse 

pinning force), as I briefly discussed here. 

• Although not discussed in these notes, noise measurements are also a good way to 
probe the periodic nature of the systems (see e.g. [61] and references therein). 

The electronic glasses are thus fascinating objects to play with, both theoretically and 
experimentally. Although we have some answers there are clearly many more unsolved 
questions than answers at the moment. Here are some possible directions to explore in 
the future in this field 

• What is the d.c. transport? To solve this question one needs first to understand 
the behavior of topological defects in a quantum crystal. Even for classical crystals 



36 



Thierry Giamarchi 



this question is open. At the moment there is only an answer for Luttinger hquids 
{d = 1) where topological defects do not play a role [20, 57]. 

• In connection with the d.c. transport, what is the out of equilibrium dynamics of 
such systems. What is the I — V characteristics. This is quite general question for 
pinned structures. For lack of space I will not discuss the question of the moving 
crystals here but one can find more details and further references in [17, 34] for 
classical crystals and in [21] for the Wigner crystal. 

• Since the system is a glass, he should show some of the properties that have been 
observed other types of glasses. In particular glasses show aging [62] , i.e. the result 
of a measurement depends on the time at which the measurement is done (the 
origin of time is fixed by the time at which the glassy phase has been formed), 
since a glass relax very slowly towards its ground state. The question of quantum 
aging is in infancy [63] and the experimental consequences for electronic glasses 
have hardly been worked out. 
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